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Abstract 



A new kind of Q-balls is found: Q-balls in a non-linear sigma model. Their main properties 
are presented together with those of their self-gravitating generalization, sigma model Q-stars. A 
simple special limit of solutions which are bound by gravity alone ("sigma stars") is also discussed 
briefly. The analysis is based on calculating the mass, global U(l) charge and binding energy for 
families of solutions parameterized by the central value of the scalar field. Two kinds (differing 
' by the potential term) of the new sigma model Q-balls and Q-stars are analyzed. They are found 

to share some characteristics while differing in other respects like their properties for weak central 
scalar fields which depend strongly on the form of the potential term. They are also compared 
(3JT)' with their "ordinary" counterparts and although similar in some respects, significant differences are 

found like the existence of an upper bound on the central scalar field. The sigma model Q-stars also 
contain non-solitonic solutions whose relation with sigma star solutions is discussed. 
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£2 ■ 1 Introduction 

■ Q-balls [TJ occur in a wide variety of (theoretical) physical contexts. They appear naturally in the 

minimal supersymmetric Standard Model [5J[2] as condensates of squarks or sleptons. The larger ones 
(Q > 10 15 ) can have a crucial cosmological significance as dark matter candidates [I] if they are stable 
\ or long living, or (since they carry baryon or lepton number) as a possible explanation for the baryon 

asymmetry in the universe [3J and the baryon to dark matter ratio [5j . Small Q-balls [6 J can be produced 
even more easily in high temperatures and may also be found as dark matter. They may also be produced 
in colliders for direct inspection of their interesting properties [5J[7]. See also e.g. Enqvist&MazumdarJS] 
and Dine&Kusenko [5] for further reviews. A large number of discussions of various other aspects of 
Q-balls exists already with different approaches: analytic [5J HU [T5J , mixed - analytic and numerical 
[13l [Ml [T5l [TBI [17] , numerical simulations [HI [19] for addressing more complicated issues like scattering 
(not yet in 3 spatial dimensions) and so on. 

All the above-mentioned Q-ball studies are based on the "original" flat space Q-balls. However, it 
is evident that for a large enough mass scale, gravitational effects become important and one needs to 
study Q-stars [SpJ dJJ [55] • The existence of Q-stars was demonstrated by Friedberg et al [53] and by 
Lynn [24] ■ Further studies revealed more features like the fact that gravity limits the size of Q-balls 
[25j or the properties of spinning Q-stars [16] . or made generalizations like Q-stars with non- minimal 
coupling to gravity [55] . 

Furthermore, unifying theories (typically in higher dimensions) lead frequently to non-linear sigma 
models, so Q-balls and Q-stars should be studied in these models as well. Although topological solitons 
in non-linear sigma models [57] have been studied for decades, it seems that non-topological solitons 
of the same models have received very little attention. This work will be devoted to a special type of 
those: Q-balls and their self-gravitating counterparts, i.e. sigma model Q-stars. 
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2 General Considerations 



Sigma model Q-balls and Q-stars are spherically symmetric solutions of the field equations derived from 
the action 

s = J d 4 x/M Qf^ixv^nv^) - um) + (2.i) 

where £ is a non- negative dimensionless function, which may be interpreted as a Weyl factor of 
a conformally-flat (two-dimensional) target space metric. A particularly simple system which will be 
studied here is the 0(3) sigma model [57] which corresponds to £ (\$\) = 1/(1 + |$| 2 /m 2 ) 2 . This is the 
conformal factor of a target space of S 2 with a diameter m. 

The function f (|$|) is a non- negative potential which will be chosen such as Q-ball solutions will 
exist as will be explained below. More notations and conventions: is the covariant derivative, the 
signature is (+, — , — , — ) and = d v T^ — dpT* v + •••. Q is Newton's constant and will denote 
the Einstein tensor. 

The field equations derived from the action (|2.ip are: 



and 



$* d£ $ dU 

£(|$|)V M V^$ + — — V A1 $V^ I $ + --— = (2.2) 
y 2|$| d\§\ ^ |$| 



- ^gG^ = = [(d^)*(d v $) + (d^)*(d^) - (V A $)*(V A $) g ^] + U(\$\)9v> (2-3) 

which is equivalent to 

^R„ v + if (|$|) [(^$)*a$) + (d v 3>T(d^)] - C/(|$|).g^ = (2.4) 

It is well-known that in flat space and a proper choice of potential which "contains any attractive 
interaction however weak" [2T], the "linear" system (£(|$|) = 1) has non-topological solitons stabilized 
by the global U(l) charge Q. The U(l) current density in the general case is given by the slightly 
modified expression 

j„ = -i£(|$|)(#*0„$ (2.5) 

The simplest way to obtain non-vanishing U(l) charge is to allow a uniform rotation in target space 
("field space"), 

$ = F{x k )e luJt (2.6) 

and indeed it can be proven [31] [5] that this must be the form of the field which minimizes the energy 
within the sector of a given Q in the "linear" theory. The generalization to the "non-linear" case is 
straightforward. The condition of finite charge leads to the boundary condition F(x ) — > at infinity. 

If we assume further spherical symmetry, F(x k ) will depend only on the radial coordinate r and the 
line element will take the usual form 

ds 2 = A 2 {r)dt 2 - B 2 (r)dr 2 ~ r 2 (d6 2 + sin 2 9dip 2 ) (2.7) 

Note that the e 4 "* factor indeed modifies the energy-momentum tensor but keeps it static. 
Einstein equations (|2 .4[) for Q-stars become 

= %-kQABt 2 ( w2£ ^ F2 - U{F)\ (2.8) 
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and they should be supplemented by the scalar field equation: 



B 2 ABr 2 A 2 J 2dF \ A 2 B 2 J dF 

The charge and mass are given by 



/>oo 

Q = 4%u drr 2 (B/A)£(F)F 2 (2.12) 
Jo 



lo 2 F 2 F' 2 \ „ \ 



M = 4. ^ At- ^(F) ^_ + — j + U(F) j (2.13) 

We note that the sigma Q-star can be viewed as a bound state of \Q\ elementary bosons (that is, it is 
stable against decay into free bosons) if M/m < \Q\. Without loss of generality we will assume oj > 
so we will have Q > as well. 

Actually, we have to solve a system of three differential equations: eq. (|2.1ip and only two of the 
three eqs (|2.8[) - (l2.10j) . By taking combinations of those three we get two first order equations (which 
may be obtained directly from the G M „ equations). 

It is more comfortable and efficient to introduce a dimensionless mass function M.(r) and use it 
instead of the metric function B(r) following the definition 

i 1 = 2GM(r) = 2M(r) 
B 2 r mr 

Note that M(r) is the accumulated mass up to radial coordinate r and the total mass of the Q-star 
is the limit M(oo) which we simply abbreviate by M where there is no danger for ambiguity. Since 
the solutions are localized, M(r) becomes essentially constant quite fast (this is one reason for using it 
instead of B(r)), and using the numerical radial end point instead of infinity is accurate enough. It is 
also simpler and more natural to revert to an angular field defined by 

|$| = mtan(6/2) (2.15) 

It is straightforward to rewrite the field equations in terms of A4(r) and 0(r) and actually to cast 
them in a dimensionless form which is ready for numerical solution. We get: 
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(2.16) 
(2.17) 
(2.18) 



where we use a dimensionless potential function it(O) and define x = mr, uj = uj/m and 7 = AirQm 2 . 

Actually, the flat space limit of this system with a conserved global charge have been studied already 
to some extent under the title "Q-lumps" [551 US], which refer to solutions which carry additional 
topological charge. This kind of topological solutions exists only in lower dimensionality or at most as 
stringlike in four dimensions. However, we will be able to find flat space (as well as self-gravitating) 
spherical solutions (in four spacetime dimensions) since we give up topological non-triviality. Therefore, 
our solutions may be simply regarded as sigma model Q-balls and Q-stars stabilized by the global charge 
alone. Within the context of the present discussion, the difference between Q-balls and Q-lumps is just 
a difference in the boundary conditions: both kinds require 0(oo) = (for finite charge) but Q-lumps 
exist for 0(0) = ir while the sigma Q-balls we find need 0(0) < ir/2. These boundary conditions are 
related to the potential functions and it turns out that Q-balls and Q-stars are easily obtained for a 
large family of potentials having a global minimum at = and another local one. We will use the 
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simple form u(0) = sin 2 (9)/8 — asin p (9)/p with p = 3, 4 so the potentials (whose local minimum is 
always at 9 = ir/2) are: 

sin 2 (6) a sin 3 (8) sin 2 (9) a sin 4 (9) 
U23(9) = , U 2i (Q) = — (2.19) 

We will see that the difference between the corresponding solutions will be analogous to those between 
the 2-3-4 and the 2-4-6 potentials of the "linear" system [T7] , 

The first term in both potentials is just a simple mass term which adds up with the u> term as in the 
"linear" system. The normalization is such that m is still the mass of the elementary free scalars. We 
choose representative values of a = 0.35 for the 2-3 potential and a — 0.4 for the 2-4 one. The potentials 
are shown for these values in figure [1] As usual, one may get the main properties of the Q-balls from 
the "effective potential" u e //(9) = u(9) — tZJ 2 sin 2 (9)/8. The parameter To will take values between 
ZZ>_ = yl — 8a/p and ZJ + = 1 and the corresponding central fields are 9,(0) = tt/2 and 9(0) — > . 

Due to the "north-south" symmetry which is left in the potential functions there exists of course 
another family of "mirror" solutions with the boundary condition 9(oo) = tt instead of the usual 
9(oo) = that we are using. The second boundary condition will satisfy accordingly 9(0) > tt/2. 
9bviously, our choice does not lead to any loss of generality of the results that will be presented here, 
but it should be kept in mind that they apply to two different families of solutions. 



3 Properties of the Solutions 

We study solutions with both potentials for the three values 7 = 0, 0.02 , 0.2. As can be guessed 
from inspection of the "effective potential", the "thin wall" and "thick wall" solutions exist for both 
potentials in a way similar to the "linear" case. 

It is obvious that there are no flat space solutions with a = 0, that is a mass term only. Gravity 
changes the situation and allows solutions which are the sigma model analogues of the boson star 
solutions, so we may call them "sigma stars" . Gne difference with respect to boson stars is that these 
sigma stars do not enjoy a scaling symmetry and do not fall on a one-parameter curve in parameter 
space. Therefore the 7 = 0.2 family of pure self-gravitating sigma model solutions which we present is 
not universal. 

The main results are shown in figures [2][6] The first, figure [2] shows the charge as a function of 
the central field 9(0) in all cases. It covers the region of maximal charge for 7 = 0.02, although the 
maximum is so narrow that its width cannot be seen it the plots. 

It may be expected that curves of the mass as a function of 9(0) are superposed in these plots, but 
they are not easily identified in this resolution. However, the plots for the binding energy per particle, 
figures [3]|4] give enough information about the masses being larger or smaller than mQ for a given 9(0). 
More insight into the situation is added by the plots of the binding energy per particle vs. charge, 
figures EM 

A new characteristic with respect to the "linear" solutions is the limiting value of the central scalar 
field 9(0) which is the same for both potentials: 9*(0) = 7r/2. Moreover, now the same limiting 9(0) 
appears also for the self-gravitating solutions and not for Q-balls only. The reason for this is the special 
form of the potential functions that we chose, which yields the minimum of the effective potential to 
be at 9(0) = tt/2 for all values of u>. We stress also that 9 = n/2 is the equator of S 2 and does not 
correspond to |$| — > 00 - see eq. (|2.15[) . Therefore, it is possible that a different choice of potential 
function will allow larger values of 9,(0). 

Another new characteristic is the existence of "mirror" solutions with different boundary conditions. 
The possibility of coexistence of solutions of both type and their interaction deserves further study. 

Generally, we find a close parallelism between the main properties of the sigma Q-balls/Q-stars and 
the analogous "ordinary" Q-balls/Q-stars. The "damped oscillations" of Q beyond its maximum which 
exist in the "linear" model [T7] exist also here but they are only partially visible in figure [5] due to the 
9(0) = 7r/2 limit. As before the solutions are stable against decay into free bosons only in a limited 
region of 9(0) values which corresponds to positive binding energy, or M/m < Q. Note however that 
charge degeneracy (two masses or more for the same Q) exists, so even in this region a higher mass state 
may decay into a lower state plus additional free bosons while conserving Q. The other solutions with 
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negative binding energies may either decay completely into free bosons, or will have a smaller stable 
boson star among their decay products. 

The sigma Q-balls with the 2-3 potential are always stable with ever growing charge. Gravity 
imposes maximal charge values: for 7 = 0.2 the maximum is at Q — 382.97 and 6(0) = 1.56921 or 
20(O)/7r = 0.99899, while for 7 = 0.02 the maximum is at Q = 32155.3. The corresponding central 
scalar field is so close to tt/2 that we give it as 1 — 29(0)/7r = 1.5543 x 10 -15 . 

The sigma star solutions with 7 = 0.2 that appear in figure [5] have a maximum at Q — 95.07 and 
O(0) = 1.54355 or 2O(0)/7r = 0.98265. These values define the region of stable solutions. As figures [3] 
and [5] show, all the others are either unbound (with negative binding energy), or may decay into the 
stable solutions while conserving particle number. 

The 2-3 Q-star curve with 7 = 0.02 in figure [2] is very similar to the 7 = Q-balls (below the 
maximum of Q — 32155.3) and it is impossible to distinguish between them. The difference shows up 
in the binding energy which is shown in figures [3^ and O 

The solutions with the 2-4 potential exhibit a more involved structure which may be described again 
according to the central field value. The small 9(0) sigma Q-balls (no gravity: 7 = 0) are large and 
unstable. The stability region starts at O(0) = 1.11515 for which Q = M/m = 23.43 and extends all 
the way to 8(0) = tt/2 with monotonically increasing charge and mass. For 7 = 0.02 the small 6(0) 
behavior changes completely. The small 6(0) solutions have positive binding energies for 6(0) < 0.14476 
(for which Q = 43.37), passing through a local maximum of charge at 6(0) = 0.06400 and Q = 50.17. 
An instability region follows for 6(0) < 1.08417. The second range of bound solutions starts at 6(0) = 
1.08417 for which Q — M/m — 21.99 and goes toward 6(0) — > tt/2. There is however a global charge 
maximum of Q = 12586.03 very close to 6(0) = tt/2. i.e. such that 1 - 26(0)/tt = 1.7686 x 10~ 13 . We 
may therefore conclude that for this 2-4 potential there are stable sigma Q-stars for any charge up to 
a maximal value of Q = 12586.03 for the parameters we chose. Some charge intervals exhibit charge 
degeneracy, so the higher mass solutions will decay. 

The solutions for 7 = 0.2 have similar behavior with two main quantitative differences: the maximal 
charge is now much smaller and has a value of Q = 200.18 and all solutions have in this case positive 
binding energy due to the stronger gravitational self-attraction. 

As far as the fiat space limit is concerned, the 2-3 solutions may be divided into two types in 
accordance with two regions of the central field interval < 6(0) < tt/2. The sigma Q-stars in most of 
this interval are similar to the flat space ones and the limit 7^0 (keeping 6(0) fixed) gives well-behaved 
sigma Q-balls. These Q-stars may be classified as solitonic since their existence does not depend upon 
gravity. The exception is the small region near (and below) 6(0) = tt/2 where there are only self- 
gravitating solutions so they are non-solitonic. The 2-4 solutions may be divided into three types: in 
the region of small 6(0) and near 6(0) = tt/2 gravity gives rise to (non-solitonic) solutions with no flat 
space limit, while in the medium 6(0) values the solutions may be considered solitonic. 

4 Summary and Outlook 

We found new Q-ball and Q-star solutions in the sigma model system. We presented the main properties 
of the solutions in flat spacetime and for two values of the gravitational strength 7 = 0.02 and 7 = 0.2. 
The corresponding characteristics are quite different for the different values of 7 as is evident from the 
plots of charge vs. central field described above, together with the analysis of the 6(0) and Q dependence 
of the binding energy. 

We found that the sigma model Q-ball and Q-star properties depend strongly on the form of the 
potential term and that the Q-star solutions split into two main types: one is self-gravitating version of 
the flat space Q-balls, while the other may be identified as non-solitonic solutions which do not have a 
flat space limit. 

There is a close parallelism with the "ordinary" Q-ball and Q-star solutions [T7] : the correspondence 
is between the sigma model with the 2-3 potential and the "linear" system with 2-3-4 potential, and 
sigma model with the 2-4 potential and the "linear" system with 2-4-6 potential. However, there are 
some important differences. The most interesting is the existence of two "mirror" families of solutions 
differing by their boundary conditions. The possibility of coexistence of solutions of both types and 
their interaction was not discussed here but deserves further investigation. 



5 



Another issue which calls for a further study is that of spinning Q-balls and Q-stars. Although 
interesting results [TBI, [T5] already exist for the "linear" system, an analogous study for the new sigma 
model solutions does not exist. A further more systematic analysis is needed in order to clarify questions 
like the relation between charge, mass and angular momentum of spinning Q-stars. 

Also needed are a better understanding of the dynamics of instability and decay processes of Q-stars 
and a study of their possible gravitational collapse. 
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Figure 1: Plots of the 2-3 potential for a — 0.35 and 2-4 potential for a = 0.4. 
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Figure 2: Plots of log(Q) vs. 9(0) for 7 = (sigma model Q-balls - dashed), 7 = 0.02 and 7 = 0.2. 
(a) 2-3 sigma model Q-stars; (b) 2-4 sigma model Q-stars. The 7 = line cannot be resolved from the 
7 = 0.02 one in the 2-3 potential. The dotted lines correspond to sigma stars with 7 = 0.2. 
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Figure 3: Plots of binding energy per particle (mQ — M)/mQ vs. 0(0) for 2-3 sigma model Q-stars with 
7 = (sigma model Q-balls - dashed), 7 = 0.02 and 7 = 0.2. (a) B.E./Q up to 0.30; (b) magnification 
of the large field region with larger B.E./Q. The dotted lines correspond to sigma stars with 7 = 0.2. 
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Figure 4: Plots of binding energy per particle (mQ — M)/mQ vs. 9(0) for 2-4 sigma model Q-stars with 
7 = (sigma model Q-balls - dashed), 7 = 0.02 and 7 = 0.2. (a) B.E./Q up to 0.15; (b) magnification 
of the large field region with larger B.E./Q. The dotted lines correspond to sigma stars with 7 = 0.2. 
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Figure 5: Plots of binding energy per particle (mQ — M)/mQ vs. log(Q) for 7 = (sigma model Q-balls 
- dashed), 7 = 0.02 and 7 = 0.2 and for boson stars with 7 = 0.2 (dotted). The insert is a magnification 
of the upper right corner which contains the 7 = 0.02 curve. 




Figure 6: Plots of binding energy per particle (mQ — M)/mQ vs. log(<2) for 2-4 sigma model Q-stars 
with 7 = (sigma model Q-balls - dashed), 7 = 0.02 and 7 = 0.2 and for boson stars with 7 = 0.2 
(dotted), (a) general view; (b) magnification of the small B.E./Q region . 
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